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A procedure is discussed for creating coherent superpositions of motional states of ion strings. The 
motional states are across the structural transition linear-zigzag, and their coherent superposition 
is achieved by means of spin-dependent forces, such that a coherent superposition of the electronic 
states of one ion evolves into an entangled state between the chain's internal and external degrees 
of freedom. It is shown that the creation of such an entangled state can be revealed by performing 
Ramsey interferometry with one ion of the chain. 
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I. INTRODUCTION 

The superposition principle is one of the fundamen- 
tal elements of quantum mechanics, most fascinating as 
it defies the everyday experience of physical reality. A 
paradigmatic example is Schrodinger's Gedankenexperi- 
ment, which pushes the quantum-mechanical predictions 
to their ultimate consequences, by envisioning an exper- 
iment that would lead to a quantum superposition in- 
volving a photon and a cat [1-3 . The "cat paradox" 
has been extensively discussed in the literature. The re- 
alization of quantum superpositions of physical objects 
having size bigger than a single atom or photon has been 
experimentally pursued in several milestone experiments, 
some of which are reported in Refs. ^4nl0 . These efforts 
ultimately aim at controlling the quantum dynamics of 
increasingly large systems, applying mostly a bottom-up 
approach. The goal is to warrant scalability of quan- 
tum mechanical operations, a basic requisite for quantum 
technological applications. 

Strings of trapped ions are one of the most promising 
platforms for quantum technologies. Among several re- 
markable achievements, we mention here that the collec- 
tive vibrations of small arrays of ions have been cooled to 
the zero-point energy [11], quantum teleportation of an 
ion's internal state has been realized p!2HT4] , entangled 
states of the internal excitations have been demonstrated 
[151 m] • First realizations of quantum simulators with 
these systems have been reported most recently [T7I42T] . 
This progress is mostly based on scaling up techniques, 
whose efficiency has been demonstrated for few particles, 
and points towards the need of identifying novel tools, 
which can allow one to control the quantum dynamics of 
mesoscopic systems. 

In this article we propose a scheme for creating 
quantum superpositions of two different structures of a 
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FIG. 1: (color online) Scheme for the creation of coherent 
superpositions of different motional states, (a) The ions are 
initially in the internal state \g) and the collective motion is in 
the ground state of a zigzag structure. A resonant laser pulse 
(green arrow) prepares the central ion in the superposition 
^ (1^) + |e)). (b) A state-dependent potential, acting only 
when the ion is in state |e), induces a conditional dynamics 
such that the ions internal and external degrees of freedom get 
entangled. When the state-dependent potential is sufficiently 
tight, the excited component will start oscillating around the 
equilibrium positions of the linear chain. 



trapped ion crystal: the linear and the zigzag order [22] . 
The superposition is accessed by driving the electronic 
transition of one ion in the chain, in a setup where an 
external field makes the trap frequency spin-dependent. 
The setup is sketched in Fig. [l] The corresponding sta- 
bility diagram is derived for three ions and a measure- 
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ment scheme, to reveal the presence of entanglement, is 
discussed, which is based on Ramsey interferometry. 

Before we start, we note that the structural transition 
from a linear to a zigzag chain has been extensively stud- 
ied in the literature [23-26 . In Ref. [24\ it was shown that 
in the thermodynamic limit this is a second order phase 
transition, which is classically described by the Landau 
model and whose control fields are the transverse trap 
frequency and the interparticle distance. If the system 
is quenched, the density of defects is well reproduced by 
the Kibble-Zurek model [27l[3Q] . Our study is a first step 
towards creating a quantum superposition of two phases 
across the quantum phase transition, while the quenching 
here is realized by coupling to a quantum degree of free- 
dom: a two-level transition of one ion in the chain. In this 
respect, our proposal extends to the ion-trap setup pre- 
vious studies which considered this idea for spin chains 
coupled to a qubit [ST-'SS . 

This article is organized as follows. In Sec. |n| we in- 
troduce the setup and the corresponding model. In Sec- 
tion [ni] schemes for preparing quantum superpositions 
of crystalline structures are proposed, while the details 
on how state-dependent crystalline structures can be cre- 
ated are reported in Sec. [TV] Measurement protocols for 
detecting the quantum superposition are discussed in 
Sec. |V| The conclusions are drawn in Sec. |VI[ while the 
Appendix provides further details on the stability dia- 
grams. 



II. THEORETICAL MODEL 

We consider N atomic ions of mass m and charge 
which are confined by an anisotropic harmonic potential 
(generated by a linear Paul trap [34] or a linearized Pen- 
ning trap [35]). We denote by vj and pj the position 
and canonical conjugate momentum of the j-th ion, with 
j = 1,...,A/'. The internal degrees of freedom of the 
ions, which are relevant to the dynamics, are the elec- 
tronic states 1^) and |e), which we assume to be stable. 
They are driven by external laser pulses, which manipu- 
late the quantum state of each ion's two-level transition 
in a controlled way. 

The total Hamiltonian H describing the dynamics of 
the ions' internal and external degrees of freedom (ex- 
cluding the laser pulses) can be decomposed into the sum 



H = i^at + ^kin + V , 



(1) 



where i^at = EjLi ^o\e)j{e\ and i^kin = EjLi Pj/(2m) 
are the internal and the kinetic energy of the ions, respec- 
tively, while V is the mechanical potential. The latter 
takes the form 
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V = Vboul + Vpot : 
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(2) 
(3) 



the unscreened Coulomb repulsion, with eo the vacuum 
permittivity. The term Vpot denotes an external, har- 
monic potential, that is assumed to be anisotropic. In 
the rest of this work we take that the confinement along 
the z-direction is so tight, that the motion along this di- 
rection can be considered to be frozen out. In the x-y 
plane, the potential depends on the internal state of the 
ions and takes the form 



N 



where 



Ve{Yj) = ]^mvl [x] + (a + Safy'^j 



(4) 

(5) 
(6) 



Here Vx is the trap frequency along the x-axis, a = fy/vx 
is the aspect ratio between transverse and axial trap fre- 
quency when the ion is in the internal state |^), and a^5a 
is the aspect ratio when the ion internal state is |e). 

The potential Q can be obtained by superimposing an 
optical potential to the ion trap. Such potential can be 
generated by a laser which is far detuned from a dipole 
transition between the state |e) and an additional auxil- 
iary electronic level, while it does not affect the dynam- 
ics when the ion is in state \g) [36, 3T. Alternatively, 
the state-dependent potential can be realized by means 
of magnetic fields, by appropriately coupling a magnetic 
dipole transition [19 . 

In this article we assume that the ions are sufficiently 
cold, such that they are localized around the equilibrium 
positions of the total potential, composed of "l^coui + "^^ot • 
The ordered structures they form depend on the the po- 
tential, and thus also on the ions' internal state. The 
potentials we will discuss support the linear and zigzag 
structures, which are illustrated on the top of Fig. [2j 

For later convenience, we introduce dimensionless vari- 
ables: the lengths are scaled by a characteristic scale 

^ = 9^/V(47^eom^J)^/^ 

which is the typical interparticle distance along the chain 
axis, while the energies by the scale 8 
dimensionless potentials read 
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where for simplicity we keep the same notation for dimen- 
sionless variables as the one we used before the rescaling. 
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III. CAT STATES OF ION CHAINS 



B. Ramsey interferometry 



In this section we describe a possible scheme for 
preparing and measuring a coherent superposition of dif- 
ferent motional states by manipulating the ions with 
lasers. We are interested in the case in which the chain is 
composed by few ions, for instance, three. We assume 
that the superimposed optical potential is sufficiently 
tight such that, when an ion in the zigzag chain is ex- 
cited, the structure becomes linear. In this regime the 
cat state is achieved by preparing a single ion in a su- 
perposition state. The parameters required in order to 
access this regime will be analysed in the next section. 



A. Creation of cat states 

The initial state of the crystal is assumed to be, un- 
less otherwise stated, given by all ions in their electronic 
ground state and the chain in the ground state of the 
zigzag structure, which we denote by |^(0)) = lO)^^- A 
laser pulse addresses one of the ions in the chain and 
drives the two-level transition \g) |e), performing the 
rotation \g) |^)cosB/2 + |e) sin 6/2. Assuming that 
the duration of the pulse is much shorter than the typical 
scales of the vibrational motion, the state of crystal after 
the pulse is |^) = |0)^^, with (for 6 = 7r/2) 



^ ( \ggg) + \geg) 



V2 



(10) 



After the pulse, the motional and internal degrees of free- 
dom get entangled by the unitary evolution governed by 
Hamiltonian Eq. ([T]), creating states of the form 



m) = ^{\999)\0).. + \9e9)\m) 



with 
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(11) 



(12) 



and He = {e\H\e) the Hamiltonian projected over the 
excited state. Here, the energy of state \g) \0)^^ is set 
equal to zero, and the state is in the reference frame of 
Hamiltonian i^at- Entanglement between internal and 
external degrees of freedom is achieved for times over 
which the overlap 



m = {m\m 



(13) 



has modulus smaller than unity. This entanglement can 
also be interpreted in terms of which-way information. 
Entanglement between internal and external degrees of 
freedom, in fact, diminuishes the visibility of the Ram- 
sey signal introducing a "distinguishability" (which-way 
information) in the interferometric path. Maximal distin- 
guishability, and correspondingly zero visibility, is found 
for / = 0. We refer the reader to Refs. ^SJIS^ for a more 
detailed discussion. 



The overlap I{t) can be measured by means of Ramsey 
interferometry, according to the scheme proposed in Ref. 
^38] and in |40| to study spin-dependent dynamics in ion 
chains. Let us assume that, after obtaining state (11) for 
a given evolution time we apply a — 7r/2 pulse, i. e. 
the inverse operation of a 7r/2 pulse. The resulting state 
reads 



1 r 
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(14)" 

and the corresponding probability to find the addressed 
ion in state |^) is given by 



Pii9) 



l + Re{J(t)} 



(15) 



Alternatively, after obtaining the state (11), we can first 



apply a phase gate on the addressed ion, namely an op- 
eration that maps the internal states according to the 
rule: \g) |^), |e) — ^ After an elapsed time t, the 
— 7r/2 pulse is applied and the probability to find the ion 
in 1^) reads 



P2{9) 



l + lm{I{t)} 



(16) 



From these two measurements the overlap I{t) is ob- 
tained. We note that this quantity is sometimes called 
"Loschmidt echo", or also "quantum fidelity" and has 
been studied extensively in other systems, for instance in 
Refs. [3I1I4TH43]. 



C. Discussion 

Some remarks are now in order. First of all, in the 
example provided here it is the ion in the middle of the 
chain which is addressed by the laser pulses. The reason 
for this choice is rather natural, as the type of excitation 
does preserve the symmetry by reflection about the cen- 
ter of the trap, which characterises both linear and zigzag 
chain. Control of the internal excitation is a necessary 
condition for the implementation of the protocol. The se- 
lective driving of the central ion could be implemented by 
a focused laser beam, as realised for instance in Ref. [13] , 
or by means of a magnetic field gradient, which tunes 
only the central ion's transition in resonance with an ex- 
ternal microwave field [44J. In this regime mechanical 
effects of the coupling to radiation can be taken to be 
very small and can be neglected [60^. The Hamiltonian 
for the pulse reads 

Hl = ^(|5)(e| - \e){g\)0mtpuise - t) (17) 

in the reference frame rotating at the laser frequency; 
9 it) is the Heaviside function, describing the step-like 
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form of the pulse with duration tpuise- Here, Q is the 
Rabi frequency, and the rotation angle of the dipole is 
0/2 l^tpuigg. 

Neglecting the evolution of the vibrational motion dur- 
ing the pulse is justified provided that the largest vibra- 
tional frequency scale, i^max, is such that i^max^puise <^ 1- 
Considering that ^t^^ise = 7r/4 for the example here con- 
sidered, the requirement can be rewritten in terms of the 
inequality i^max ^ ^• 

We note that the scheme here proposed is based on cre- 
ating a local defect when the central ion is in the excited 
state: this defect will induce another ordering (e.g., when 
starting from a zigzag array, the ion excitation induces 
a transition to a linear chain) provided that the length 
of the string is shorter than the correlation length [24] . 
However, as the number of ions is scaled up, a linear su- 
perposition like the state in Eq. ( 10) will present a central 



deformation of finite size. This problem could be solved 
by creating a Greenberger-Horne-Zeilinger (GHZ) state, 
namely. 



1 
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ilggg) + |eee)) , 



(18) 



which warrants that all ions feel the same potential. This 
kind of state has been realized in [45*, "46^ [GT. The prepa- 
ration of this class of states becomes increasingly difficult 
for larger ion chains [45l |46], and the entanglement of 
GHZ-states is in general not robust [47] . 

To conclude this section, we mention that an alterna- 
tive procedure for generating a cat state could be imple- 
mented by (i) first preparing the superposition state ([lO| 
for the internal degrees of freedom, and then (ii) slowly 
switching on the additional state-dependent potential, 
such that the vibrational state follows adiabatically the 
change in the (state-dependent) vibrational Hamiltonian. 
In this form, it would be possible to prepare a superposi- 
tion of two different structures in their respective ground 
state. This proposal poses several challenges and it will 
be elaborated in more detail in a subsequent publication. 



IV. STATE-DEPENDENT CRYSTALLINE 
STRUCTURES 

In this section we shall analyse the equilibrium configu- 
ration and normal modes of the ion chain in the harmonic 
trapping potential. The stationary equilibrium positions 



.(0) 



satisfy the equations |^(r^-^^) = 0, and they are 
stable when the eigenvalues of the Hessian matrix of V 
evaluated at r^-^^ are larger than zero. In this case the 
eigenvalues give the frequencies of the normal modes and 
the corresponding eigenvectors are the eigenmodes [48] . 
Clearly, in our case there will be different sets of solutions 
depending on the ions' internal state. 

In the following we will consider that the ion excitation 
couples the ground states of the classical phase transition 
|24|. We will hence discard the regime, in which one 



would observe the disordered phase, predicted for the 
quantum phase transition, and which is observed when 
the trap frequency is close but below the critical value, 
yet quantum ffuctuations destroy the zigzag order [26^, 
^49^. This would correspond to a narrow interval of values, 
that has been characterised in detail in Ref. [49 . In 
Ref. ^5 a characterisation of the parameters required 
for small chains of ions can be found. 



A. Spatial organization of the ions in harmonic 
potentials 

We first examine the equilibrium structures for three 
trapped ions in the case in which the internal state of all 
particles is \g) and the confinement in the ^/-direction is 
tighter than that along x, i.e., a > 1. The following dis- 
cussion reviews results which have been reported, for in- 
stance, in [24 , 50 , 51 . For this case there are two different 
possible types of solutions for the equilibrium positions: 
the linear chain, with all ions aligned along the x-axis, 
and a planar structure, where the ions form the extremes 
of an isosceles triangle whose symmetry axis is aligned 
along the ?/-axis. We refer to this planar structure as 
zigzag configuration. More specifically, the equilibrium 
positions of the linear chain are 
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0, 



(19) 
(20) 



T his co nfiguration is stable when a > ac, with ac = 
^"12/5 ^ 1.5492. In this case, the axial normal mode 
frequencies are {1, a/3, y^29/5} and the transverse ones 
are {\/a'^ — Vc?^^^, a}. 

For a < ac (and a > 1) the chain is in a zigzag con- 
figuration, with equilibrium positions 
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0, 
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where 



a" 

3 \2/3 



-1/3 



.1/2 



(21) 
(22) 

(23) 
(24) 



The analytic expressions for the normal mode frequencies 
are rather cumbersome and are reported, for instance, in 
Ref. [24 for the case of N ions in a ring. Figure [5] displays 
the frequencies of the normal modes as a function of the 
aspect ratio a. 



B. Spatial organization in state- dependent 
harmonic potentials 

We now extend the previous considerations to the case 
in which the internal state of the three ions is not the 
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FIG. 2: Normal mode frequencies for three ions confined in 
a harmonic potential, in units of the trap frequency z/^, as a 
function of the aspect ratio a — Vyjvx- The dashed vertical 
line indicates the critical value cxc ~ 1.5492, separating the 
zigzag [a < ac) from the linear array {a > ac). 



same, so that the ions experience different trapping po- 
tentials. We shall consider only the case in which one of 
the three ions is excited. In particular, we will discuss the 
case in which the excited ion is in the middle of a linear 
chain. The configurations found for the case, in which 
the excited ion is at one edge of the chain, are reported 
for completeness in the Appendix. 

When the ions form a linear array, the equilibrium con- 
figurations corresponding to state \geg) are the same as 
in the homogeneous case, and are thus given by Eqs 



(19 



20). In this case, in fact, the central ion is located in the 
center of the trap where the external potential is zero. 
For this case the eigenfrequencies read 

{l,V^, ^2975, - 1, 

^/a^ + aSa - (12 - bSa'^ - p)/10, 

+ aSa - (12 - 5fo2 + p)/10 | , (25) 



with p = V'128+ [5(5a(2a + Ja) -4]2. By imposing that 
the eigenfrequencies are real and positive, one finds that 
the linear chain is stable for 5a > Sac^ with 



Sac = 




l a. 



(26) 



where a > 1. The corresponding curve is reported in the 
stability diagram of Fig. [S] Here, one observes that for 
most values of Sa < Sac the chain is in a zigzag structure 



along the x-axis (i.e. symmetric about the y-axis), whose 
equilibrium positions read 
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and 



R 



{a + Sa) 



2 • 



(27) 
(28) 

(29) 
(30) 

(31) 



The homogeneous case is recovered in the formula by 
setting Sa = 0: in this limit, R = 1 and Eqs. ([29|, ([30| 
coincide with Eqs. (23), (24), respectively. 

The region in the upper left corner of Fig. |3] corre- 
sponds to the appearance of zigzag configurations along 
the ^/-direction. The corresponding stability boundary 
is evaluated by imposing that the Hessian matrix, giv- 
ing the normal modes, has vanishing determinant, which 
leads to the equation: 



^ _ (y_y 

2R^l \dJ 



{3xy) 



(32) 



where x, y and D correspond to the equilibrium positions 
in this configuration and are given by: 



y = 
D 



1/2 



(9^2 + y 



1/3 

:2^1/2 _ „i/3 



(33) 
(34) 
(35) 



We note that, as opposed to the case when all ions expe- 
rience the same potential, more than one configuration 
can be stable in the same region. 



V. LOSCHMIDT ECHO 

A clear demonstration of the creation of cat states 
would require quantum tomography of the crystal state 
j8[ [521154 . Nevertheless, signatures of the quantum su- 
perposition can be obtained by means of Ramsey interfer- 
ometry. This scheme allows one, as we have shown, to de- 
termine the overlap I{t) as a function of the elapsed time 
t between the two laser pulses. In this section we analyse 
the behaviour of the absolute value of the overlap \I{t)\ 
for different values of the axial and spin-dependent poten- 
tial. This is performed by scanning through various re- 
gions of the phase diagram of the linear-zigzag structural 
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FIG. 3: (color online) Stable configurations for a three-ion 
chain as a function of the aspect ratio a and variation 5a (here 
due to a state-dependent potential acting on the central ion). 
The labels "LIN X" , "ZZ X" and "ZZ Y" correspond respec- 
tively to a linear chain along the x-axis, a zigzag chain along 
the X-axis and a zigzag along the ^/-axis. The solid lines are 
the boundaries given in Eqs. (26) and (32). The other bound- 



ary was determined numerically (the wiggled boundary is here 
due to the chosen numerical resolution) . The variation of the 
state-dependent potential corresponds to a vertical displace- 
ment in the plot. 



instability for the case of three ions. For this calculation, 
the potentials were approximated by their quadratic ex- 
pansions about the corresponding equilibrium positions, 
and the vibrations have been quantized according to the 
procedure reported for instance in Ref. [50 . 

In the following discussion we will assume that the 
chain is initially either in the zigzag or in the linear con- 
figuration, depending on the value of the aspect ratio a, 
and that the ions are prepared in state (11) at t = 0, 



as described in Sec. |III[ The strength of the harmonic 
potential trapping the central ion is given by the fre- 



quency v[ 



y 



SUy^ with SUy 



Sa I 



We first focus 



on the signal for different values of a, fixing the quench 
Sa. This corresponds to analysing the Ramsey signal 
by keeping constant the frequency shift dvy (due to the 
spin-dependent potential) but varying the potential at 
the electrodes. 

Figure |4] displays \I{t)\ for different values of a when 
the excitation of the central ion gives rise to a transi- 
tion (a) between two zigzag structures (where the second 
one is more tightly confined), (b) from a zigzag to a lin- 
ear string, and (c) between two linear strings (where the 
second one is more tightly confined). We note that the 
observed time-dependent behaviour is determined by the 
modes of the crystalline structure that is generated by 
exciting the central ion. In fact, the initial structure, 
which corresponds to the structure when all ions are in 
state 1^), is the ground state of the corresponding array. 
In (a) revivals of the overlap signal are observed. These 
revivals become more distant in time the closer the zigzag 
array is to the mechanical instability, where the time in- 
tervals scale with the period of the lowest eigenfrequency 
of the tighter zigzag chain. In (b), where the excitation 
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FIG. 4: Overlap signal |/(t)| as a function of the elapsed 
time t (in /is) for a Ramsey setup in which the central ion 
of a three-ion chain is excited. The curves are evaluated for 
^Be^ ions for the parameters zy^ = 27r x 500 kHz, i/y = avx, 



and Uy = ly^ 



■y -\- 27T X 10 kHz (corresponding to Sa = 0.02). 
The values of a (and thus Uy) are given in the boxes. The 
equilibrium structures for ground and excited states are (a) 
zigzag-arrays with different transverse displacements (here 
Uy = 27T X {725,745,765} kHz); (b) a zigzag- and a linear 
array, respectively (uy = 27r x {772.5,773.0,773.5} kHz); (c) 
linear arrays with different transverse confinements (i^y — 
2ti X {775.0,776.0,780.0} kHz). 
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corresponds to a quench across the transition connect- 
ing zigzag and hnear chain, one stih observes revivals of 
the overlap, whereby the signal exhibits a more complex, 
quasi-periodic structure. When the transition is instead 
between two linear arrays with different transverse con- 
finement (c), the overlap exhibits a periodic modulation 
whose maximal amplitude is unity and at twice the fre- 
quency of the lowest frequency eigenmode: the zigzag 
mode of the linear chain. This modulation is a manifes- 
tation of squeezing induced by the sudden change of the 
trap potential. 

This behaviour can be further characterised using the 
Fourier spectra of the overlap signal. These are shown 
in Figs. [5] for each of the three cases displayed in Figs.|4] 
A first look shows that the spectral components of the 
transition between the zigzag structures are dominated 
by a comb of frequencies. Fig. |4ja), with a background 
signal. The comb is at frequencies that are multiple of 
the lowest frequency of the tight zigzag array. This struc- 
ture splits up into more components when the transition 
couples two different structures. Fig. [4]^b). Finally, when 
the transition couples two linear structures, Fig.[4]^c), the 
spectrum is much more regular and composed by only 
a few frequency components, dominated by the lowest 
one. Further analysis shows that the dominant contribu- 
tion in Fig. |4]^a) and (c) is from the zigzag mode, which 
is also the soft mode of the structural phase transition. 
This mode seems to play a relevant role also when the 
transition couples two structures across the structural in- 
stability. Fig. Qb). The revivals are associated with the 
oscillations of this mode, which is mostly excited by the 
quench. 

We note that the appearance of the revival signal is 
a signature of the quantum superposition: this signal 
would be absent if there is no quantum coherence be- 
tween the two crystalline states. We then focus on its 
features, and characterise them in detail considering elec- 
tronic excitations that couple a crystal with an equilib- 
rium zigzag structure to a crystal with a linear array. 

The figures discussed in the following refer to the sit- 
uations in which ground and excited states correspond 
to equilibrium structures across the mechanical insta- 
bility, namely, when the initial state is a zigzag array 
and the excited state is a linear chain. Figure [6] shows 
the overlap signal for different frequencies of the spin- 
dependent potential, corresponding to different values of 
the parameter 5a (larger absolute values correspond to 
larger variations). The distance between the revival sig- 
nals decreases as da increases, namely, the excited state 
is deeper in the linear phase. Moreover, for Sa = 0.1 one 
observes the onset of a periodic signal. The form of the 
signal at da = 0.02 is further characterised in Fig. [7| as 
a function of the axial trap frequency, showing a slight 
increase of the revival signals as is increased, and cor- 
respondingly as the ions are more tightly confined in the 
x-direction. 

Finally, the form of the revivals is reported in Fig. [S] as 
a function of the ion species, that is, of the mass deter- 
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FIG. 5: Fourier transform of the overlap signal |/(t) | in Fig. [4] 
Here, (a) corresponds to the curve at a = 1.49, (b) to a = 
1.547, and (c) to a = 1.55. 



mining the spread of the motional wave packet. In this 
case, for lighter ions, i.e., for larger quantum fluctuations, 
one observes a slight increase in the revival signal. 

The dynamics we considered so far is unitary: decay 
of the overlap signal is solely due to dephasing in the dy- 
namics of internal and external degrees of freedom. De- 
coherence and noise sources have been neglected. These 
would affect both internal and external degrees of free- 
dom. They are expected to introduce a damping factor in 
the overlap signal, setting an upper bound for partial re- 
vivals decreasing with time, such that a measured revival 
in an experiment in turn would give an estimated bound 
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FIG. 6: Same as Fig. [I] but for a = 1.5470 [vy = 27r x 
773.5 kHz) and 5a = {0.01,0.04,0.10} (corresponding to 
dfy = 27r X {5,20,50} kHz). In this case the equihbrium 
structures for ground and excited states are a zigzag and a 
hnear array, respectively. 
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FIG. 7: Same as Fig. [6] but for 5a = 0.02 and differ- 
ent axial trap frequencies Vx (the values are reported in 
the box). The value for a corresponds to z/^ = 27r x 
{154.7, 309.4, 773.5, 1547} kHz, the frequency shift due to the 
optical potential is dvy — 2ti x {2, 4, 10, 20} kHz (ordering 
corresponds to the values of Vx reported in the legenda in the 
box, from bottom to top). The color-coded abscissa indicate 
the time (in /xs) for the curves of the corresponding color. 



on decoherence effects. In order to observe the curves re- 
ported in Fig. [4j one needs coherence times longer than 
100z^~^ for both internal and external degrees of freedom. 
Taking — 27rx500 kHz, the state would need to remain 
coherent for times longer than 30 /is. This assumption 
seems reasonable given the coherence times reported for 
instance in P |55l [5^ . The requirements of our proposal 
in terms of optical power are very moderate: for the val- 
ues corresponding to Fig. [4j the optical potential should 



FIG. 8: Same as Fig. [6] but for 5a = 0.02 but for different 
ion species. 



have an associated motional frequency of the order of 
27r x 100 kHz, well within present reach. Finally, we 
note that the ability to combine optical potentials and 
ion traps has already been demonstrated [37l |57] . 



VI. CONCLUSIONS 

In this paper we have studied the dynamics of the 
structural phase transition from the linear to the zigzag 
configuration in a chain of trapped ions, with the aim 
of creating a coherent superposition of different meso- 
scopic crystalline structures. Among several possibili- 
ties to achieve that goal, we have chosen to exploit spin- 
dependent trapping potentials, which influence each ion's 
motional state as a function of its internal electronic 
state. This allows in particular to create an entangled 
state from coherent superpositions of one degree of free- 
dom. 

Specifically, our protocol starts with an ion chain 
at rest under external confinement in a regime near 
the chain's structural phase transition. Excitation of 
one particular ion from the chain into a superposition 
of metastable electronic states affects, via the state- 
dependence of the external potential, this ion's motion, 
entangling its vibrational and internal state. As a result 
of Coulomb repulsion, the other ions in turn rearrange 
their motion according to this ion's new state, thereby 
becoming entangled with it as well. 

This realizes, under appropriate conditions described 
in the text, a mesoscopic superposition (a cat state) of 
different crystalline structures, something that to our 
knowledge has not been observed or proposed before [62j . 
Signatures of the generated cat state can be obtained by 
means of Ramsey interferometry, according to the scheme 
discussed in [38 . Our analysis shows that these dynam- 
ics could in principle be observed in current experiments, 
using crystals of three ions. 

Future work will focus on the properties of the Ram- 
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sey signal measured for longer chains when the ions are 
in the disordered phase of the quantum phase transition 
[26l [49] . Moreover, the efficiency of the protocol will 
be characterised as a function of the temperature of the 
chain, including also heating and decoherence effects sys- 
tematically in the theoretical model. The control tools 
proposed here open the way to the application of optimal 
control techniques in the spirit of the proposals in Refs. 
[58l [59] , for the purpose of creating robust mesoscopic 
quantum states on demand. 
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Appendix A: Equilibrium configurations when one 
outer ion is excited 



We now determine the stability diagram when one 
of the outer ions is excited and experiences a state- 
dependent potential. We distinguish between three dif- 
ferent cases: 1) a linear chain configuration along the 
x-direction, which we denote by LIN X*; 2) a linear 
chain configuration along the y-axis, which we denote by 
LIN Y*; and 3) a triangular configuration which we label 
by TRIA* (here we choose not to distinguish between 
orientation along x or y). The asterisk always denotes 
the structures where one of the outer ions is excited. As 
in the previous cases, we shall focus on the parameter 
region a > 1, for which a linear chain along the y-ax.is is 
never stable. 

For the LIN X* we obtain again the same equilibrium 
positions as in Eq. (TgjpQ). It is possible to determine 



the stability boundary analytically by solving equation 

(65/8)^2 - 9 



2R 
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5(5a4-25aV2 + 4) ' 



(Al) 



where R is defined as in Eq. (31). 



The triangular structure TRIA* does not exhibit the 
symmetry properties of the zigzag structure. The cor- 
responding equilibrium positions, as well as the stability 
boundary shown in Fig. |9| were calculated numerically 
using a Metropolis algorithm with random initial posi- 
tions followed by a constraint minimization setting the 
asymmetry as a constraint. However, as the structures 
transform smoothly, one has to allow a certain tolerance 
on the constraint. On the other hand, the algorithm 
1 



6a ^ 




-0.5 



FIG. 9: (color online) Same as Fig. [sjwhen one outer ion 
experiences the state-dependent potential. The labels "LIN 
X*" and "TRIA*" correspond respectively to a linear chain 
along the x-axis, and a triangular structure. The solid lines 
correspond to the stability boundaries in Eqs. (32) and (Al), 
the others were calculated numerically. The white regions 
are such that there is no stable structure with one outer ion 
excited. 



might run into a minimum corresponding to another 
structure satisfying accidentally the constraints. We be- 
lieve this is the reason why scattered singular solutions 
are observed in the stability diagram for the triangular 
configuration. 



[1] E. Schrodinger, Naturwissenschaften 23, 807, 823, and 
844. (1935), English translation by J. D. Trimmer, Proc. 
Am. Phys. Soc. 124, 3235 (1980). 

[2] W.H. Zurek, Phys. Today 44, 36 (1991). 

[3] J. A. Bergou and B.-G. Englert, J. Mod. Opt. 45, 701 
(1998). 

[4] A. Auffeves, P. Maioli, T. Meunier, S. Gleyzes, G. 

Nogues, M. Brune, J. M. Raimond, and S. Haroche, Phys. 

Rev. Lett. 91, 230405 (2003). 
[5] S. Gerlich, S. Eibenberger, M. Tomandl, S. Nimmrichter, 

K. Hornberger, P. J. Fagan, J. Tiixen, M. Mayor, and 

M. Arndt, Nature Communications 2, 263 (2011). 



[6] W. Wieczorek, R. Krischek, N. Kiesel, P. Michelberger, 
G. Toth, and H. Weinfurter, Phys. Rev. Lett. 103, 
020504 (2009). 

[7] R. Prevedel, G. Cronenberg, M. S. Tame, M. Paternostro, 

P. Walther, M. S. Kim, and A. Zeilinger, Phys. Rev. Lett. 

103, 020503 (2009). 
[8] D. Leibfried, R. Blatt, C. Monroe, and D. Wineland, Rev. 

Mod. Phys. 75, 281 (2003). 
[9] H. Haffner, C. F. Roos, and R. Blatt, Phys. Rep. 469, 

155 (2008). 

[10] A. D. O'Connell, M. Hofheinz, M. Ansmann, R. C. Bial- 
czak, M. Lenander, E. Lucero, M. Neeley, D. Sank, 



10 



H. Wang, M. Weider, J. Wenner, J. M. Martinis, and 

A. N. Cleland, Nature 464, 697 (2010). 
[11] J. D. Jost, J. P. Home, J. M. Amini, D. Hanneke, R. Oz- 

eri, C. Langer, J. J. Bollinger, D. Leibfried, and D. J. 

Wineland, Nature (London) 459, 683 (2009). 
[12] S. Olmschenk, D. N. Matsukevich, P. Maunz, D. Hayes, 

L.-M. Duan, and C. Monroe, Science 323, 486 (2009). 
[13] M. Riebe, H. Haffner, C. F. Roos, W. Hansel, J. Benhelm, 

G. P. T. Lancaster, T. W. Karber, C. Becher, F. Schmidt- 
Kaler, D. F. V. James, and R. Blatt, Nature 429, 734 
(2004). 

[14] M. D. Barrett, J. Chiaverini, T. Schaetz, J. Britton, 
W. M. Itano, J. D. Jost, E. Knill, C. Langer, D. Leibfried, 
R. Ozeri, and D. J. Wineland, Nature 429, 737 (2004). 

[15] D. Leibfried, E. Knill, S. Seidelin, J. Britton, 
R. B. Blakestad, J. Chiaverini, D. B. Hume, W. M. Itano, 
J. D. Jost, C. Langer, R. Ozeri, R. Reichle, and 
D. J. Wineland, Nature 438, 639 (2006). 

[16] T. Monz, P. Schindler, J. T. Barreiro, M. Chwalla, 
D. Nigg, W. A. Coish, M. Harlander, W. Hansel, M. Hen- 
nrich, and R. Blatt, Phys. Rev. Lett. 106, 130506 (2011). 

[17] R. Islam, E. E. Edwards, K. Kim, S. Korenblit, C. Noh, 

H. Carmichael, G.-D. Lin, L.-M. Duan, C.-C. Joseph 
Wang, J. K. Freericks, et al.. Nature Comm. 2 (2011). 

[18] B. P. Lanyon, C. Hempel, D. Nigg, M. Miiller, R. Ger- 
ritsma, F. Zahringer, P. Schindler, J. T. Barreiro, 
M. Rambach, G. Kirchmair, et al., ArXiv e-prints (2011), 
1109.1512. 

[19] M. Johanning, A. F. Varon, and C. Wunderlich, J. Phys. 

B 42, 154009 (2009). 
[20] C. Schneider, D. Porras, and T. Schaetz, ArXiv e-print 

(2011), 1106.2597. 
[21] J. T. Barreiro, M. Muller, P. Schindler, D. Nigg, T. Monz, 

M. Chwalla, M. Hennrich, C. F. Roos, P. Zoller, and 

R. Blatt, Nature 470, 486 (2011). 
[22] K. S. Birkl, G. and H. Walther, Nature 357, 310 (1992). 
[23] G. Morigi and S. Fishman, Phys. Rev. Lett. 93, 170602 

(2004). 

[24] S. Fishman, G. De Chiara, T. Calarco, and G. Morigi, 

Phys. Rev. B 77, 064111 (2008). 
[25] A. Retzker, R. C. Thompson, D. M. Segal, and M. B. 

Plenio, Phys. Rev. Lett. 101, 260504 (2008). 
[26] E. Shimshoni, G. Morigi, and S. Fishman, Phys. Rev. 

Lett. 106, 010401 (2011). 
[27] T. W. B. Kibble, J. Phys. A 9, 1387 (1976). 
[28] W. H. Zurek, Nature (London) 317, 505 (1985). 
[29] A. del Campo, G. De Chiara, G. Morigi, M. B. Plenio, 

and A. Retzker, Phys. Rev. Lett. 105, 075701 (2010). 
[30] G. De Chiara, A. del Campo, G. Morigi, M. B. Plenio, 

and A. Retzker, New. J. Phys 12, 115003 (2010). 
[31] H. T. Quan, Z. Song, X. F. Liu, P. Zanardi, and C. P. 

Sun, Phys. Rev. Lett. 96, 140604 (2006). 
[32] J. Zhang, X. Peng, N. Rajendran, and D. Suter, Phys. 

Rev. Lett. 100, 100501 (2008). 
[33] C. Cormick and J. P. Paz, Phys. Rev. A 77, 022317 

(2008). 

[34] P. K. Ghosh, Ion traps (Clarendon Press, 1995). 

[35] H. F. PoweU, D. M. Segal, and R. C. Thompson, Phys. 

Rev. Lett. 89, 093003 (2002). 
[36] R. Grimm, M. Weidemiiller, and Y. B. Ovchinnikov, Adv. 

At. Mol. Opt. Phys. 42, 95 (2000). 
[37] C. Schneider, M. Enderlein, T. Huber, and T. Schaetz, 

Nature Phot. 4, 772 (2010). 



[38] G. De Chiara, T. Calarco, S. Fishman, and G. Morigi, 

Phys. Rev. A 78, 043414 (2008). 
[39] B.-G. Englert, Phys. Rev. Lett. 77, 2154 (1996). 
[40] T. Pruttivarasin, M. Ramm, I. Talukdar, A. Kreuter, and 

H. Haeffner, New J. Phys. 13, 075012 (2011). 

[41] P. Jacquod, P. G. Silvestrov, and C. W. J. Beenakker, 

Phys. Rev. E 64, 055203 (2001). 
[42] F. M. Cucchietti, D. A. R. Dalvit, J. P. Paz, and W. H. 

Zurek, Phys. Rev. Lett. 91, 210403 (2003). 
[43] T. Gorin, T. Prosen, T. H. Sehgman, and M. Znidaric, 

Phys. Rep. 435, 33 (2006). 
[44] C. Balzer, A. Braun, T. Hannemann, C. Paape, M. Et- 

tler, W. Neuhauser, and C. Wunderlich, Phys. Rev. A 

73, 041407 (2006). 
[45] C. F. Roos, M. Riebe, H. Haffner, W. Hansel, J. Benhelm, 

G. P. T. Lancaster, C. Becher, F. Schmidt-Kaler, and 

R. Blatt, Science 304, 1478 (2004). 
[46] D. Leibfried, E. KniU, S. Seidelin, J. Britton, R. B. 

Blakestad, J. Chiaverini, D. B. Hume, W. M. Itano, J. D. 

Jost, C. Langer, et al.. Nature 438, 639 (2005). 
[47] M. Hein, W. Diir, and H.-J. Briegel, Phys. Rev. A 71, 

032350 (2005). 

[48] H. Goldstein, Classical Mechanics (Addison- Wesley 
Press, 1950). 

[49] E. Shimshoni, G. Morigi, and S. Fishman, Phys. Rev. A 

83, 032308 (2011). 

[50] G. Morigi and S. Fishman, Phys. Rev. E 70, 066141 
(2004). 

[51] D. F. V. James, App. Phys. B 66, 181 (1998). 

[52] S. Deleglise, I. Dotsenko, C. Sayrin, J. Bernu, M. Brune, 

J.-M. Raimond, and S. Haroche, Nature (London) 455, 

510 (2008). 

[53] L. G. Lutterbach and L. Davidovich, Phys. Rev. Lett. 

78, 2547 (1997). 

[54] J. F. Poyatos, R. Walser, J. I. Cirac, P. Zoller, and 
R. Blatt, Phys. Rev. A 53, R1966 (1996). 

[55] L. Deslauriers, S. Olmschenk, D. Stick, W. K. Hensinger, 
J. Sterk, and C. Monroe, Phys. Rev. Lett. 97, 103007 
(2006). 

[56] N. Timoney, I. Baumgart, M. Johanning, A. F. Varon, 
M. B. Plenio, A. Retzker, and C. Wunderlich, Nature 
476, 185188 (2011). 

[57] H. Katori, S. Schhpf, and H. Walther, Phys. Rev. Lett. 

79, 2221 (1997). 

[58] P. Doria, T. Calarco, and S. Montangero, Phys. Rev. 

Lett. 106, 190501 (2011). 
[59] T. Caneva, T. Calarco, and S. Montangero, Phys. Rev. A 

84, 22326 (2011). 

[60] Neglecting the mechanical effects of the laser pulse sim- 
plifies the theoretical treatment but does not limit the 
feasibility of the protocol. 

[61] It should be noticed that its preparation usually makes 
use of collective excitations of the ion chain, such that the 
free evolution of the vibrational motion during the prepa- 
ration pulses cannot be neglected. Nevertheless, this just 
implies that the initial state of the chain is not the ground 
state of the motion but still a pure state 

[62] After completion of this work we became aware of 
work along related lines carried out in: W. Li and 

I. Lesanovsky, ArXiv e-print (2011), 1108.3591, where 
the creation of a three-ion state-dependent structure was 
proposed by means of Rydberg excitation of the central 
ion. 



